Chapter 3 
Quantum Models 



In classical statistical mechanics we did not know the exact microscopic state of 
the system (location in the phase space). We made use of the ergodic hypothesis 
and replaced the time average of a physical quantity with an ensemble average, 
i.e., an average over many equivalent systems. 

In quantum statistical mechanics the situation is similar. 



3.1 Quantum Statistics 

In this Chapter we will first introduce the density matrix and then make an intro- 
duction to the Bose condensation. 

In quantum mechanics, the energy levels are discrete and the state of the sys- 
tem is described by a wavefunction. Let us first assume that we know the dynami- 
cal state (wavefunction) \ip) of the whole system and the complete set of stationary 
states \n) which are the eigenstates of the full Hamiltonian of the compound sys- 
tem, 

U\n) = E n \n). (3.1) 

(We assume that the states \n) are normalized.) n denotes the quantum states (e.g., 
momenta) of all particles in the system. E n is the total energy of all particles of 
the system in the compound state n. Any dynamical state of the system can be 
expressed as a linear combination of these eigenstates, 

\m)=Y, c n(t)\n), (3.2) 

n 
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The coefficients c n (t) define a point in the (Hilbert) space of wavefunctions \n). 
The (time-dependent) wavefunction \tp(t)) which obeys the Schrodinger equation 

ih^- = H\1>) (3.3) 

is thus the quantum-mechanical analogue of the point in phase space of classical 
statistical mechanics. The expectation value of a physical quantity A (which does 
not necessarily commute with H) is: 

{A{t)) = (m\A\m) = £<c m (n|A|m> 

m,n 

(3.4) 

The last expression holds for any complete set of states \n), also if they are not 
the eigenstates of the full Hamiltonian. 



3.1.1 The Density Matrix 

As the number of particles in the system increases, the separation between the 
energy levels decreases rapidly, and, in the thermodynamic limit, it becomes ex- 
tremely high. On the other hand, the energy levels are never completely sharp. 
They are broadened for many reasons, including the uncertainty principle. In a 
macroscopic system (with extremely dense energy levels), the levels will thus al- 
ways overlap. Consequently, a macroscopic system will never be in a strictly sta- 
tionary quantum state, it will always be in a dynamical state, in a time-dependent 
mixture of stationary states. 

It is neither feasible nor possible to find a complete description of such a sys- 
tem. Like in classical statistical mechanics, we first assume that the system is 
ergodic, we replace the time average by an ensemble average over many quantum 
systems at the same instance. We talk about probabilities of finding the system 
in the states \n). The dynamical state is no longer a unique quantum state but 
is a statistical mixture of quantum states. We introduce the density matrix as the 
ensemble average of c* m c n : 

p = \m)p m;n {n\ 

p m ,n = (mlpln) = -Y,Cn*c l m . (3.5) 
" 1=1 

The sum over i runs over £1 configurations of the system, all being part of 
the same ensemble, and c l m = (m\ip l ). We do an ensemble average because our 
information on the system is not complete. 
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The density matrix is a statistical operator with the following important prop- 
erties: 

• The density matrix completely defines an ensemble of quantum systems 
and it carries all the information that is available for a quantum statistical 
ensemble. 

• Its diagonal elements p„ jn tell the probability that the system is in the state 
\n) whereas its non-diagonal elements p m n tell the probability of a transi- 
tion from the state \n) to the state \m). For stationary states, thus, p has to 
be diagonal, it commutes with the Hamiltonian. 

• The density matrix is normalized, 

Tr(p) = 1. (3.6) 

The trace runs over any complete set of states \n) and is independent of the 
choice of the basis set. 

• The (ensemble and quantum-mechanical) average value (A) of any physical 
quantity A is equal to: 

(A) = (ij\A\ij) = £ Pm,n(n\ A \™) = Tr(pA) = Tr(Ap). (3.7) 

n,m 

The trace in quantum statistical mechanics plays the role of integration over 
the whole available phase space in classical statistical mechanics. 

• The dynamics of the system is completely described by: 

ih^ = [H,p}. (3.8) 
This is the "Schrodinger equation" of the density matrix. 

3.1.2 Ensembles in Quantum Statistical Mechanics 

Now, we shall specify the density matrix and its relation to thermodynamic quan- 
tities. 
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Quantum Microcanonical ensemble. 

The energy of systems belonging to a microcanonical ensemble is fixed, 
therefore it is convenient to use the energy eigenstates as a basis. Then the 
density matrix is diagonal. According to the postulate of equal weights, 
all the states \n) with the energy E n between E and E + AE are equally 
probable and let there be AVI such states. For a quantum system, Ail is an 
integer number. The probability for the system being in one of these states 
is p = 1/Afi and the entropy is equal to: 



S = -k B (\np). 



(3.9) 



In general, we must use the density matrix instead of simple probabilities p 
and the entropy of a quantum system is: 



S = —k B (\np) = —ks Tr(plnp) 



(3.10) 



Again, the average ( ) means ensemble and quantum average - both aver- 
ages cannot be separated. 

Canonical ensemble. In the energy representation p is diagonal and we can 
use the same arguments as in the case of classical canonical distributions. 
The matrix elements are equal to: 

1 



Pr, 



-PE„ 



where the normalization constant is (again) the partition function 



-PE„ 



(3.11) 



(3.12) 



In an arbitrary basis, the density matrix can be written more generally as 




with 



Tre 



-PH 



(3.13) 



(p and H are operators!). Thus, for a canonical ensemble, the average (ther- 
mal and quantum- mechanical) of an operator A is: 



1 



(A) = Tr (pA) = - Tr (e"^ A) 



(3.14) 
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In particular, the internal energy is 



E=(H) = l-Tx (He-*™) = - A i nT r e"^ 
/' ^ up 



The free energy is 



E=-^- j \nZ. (3.15) 
op 



-k B T\n Z. (3.16) 



All the thermodynamic relations are the same as before (Tables I - III are 
still valid), the only difference is that now we have to calculate the partition 
function as the trace of the density matrix. 

The entropy in a canonical ensemble, S, is the ensemble average over all 
(energy) states of the microcanonical entropy S: 

S = -(k B \n p) = -fc B Tr(plnp). (3.17) 

• Grand canonical ensemble. In the grand canonical ensemble the operator p 
operates on a generalized Hilbert space which is the direct sum of all Hilbert 
spaces with fixed number of particles. The density matrix is: 



where the grand partition function is 

~(T,V,p) =Tre-^ n ^ N \ 

Notice that now iV is an operator in the generalized Hilbert space. As in the 
classical case, S is equal to the canonical partition function Z, weighted by 
the fugacity z = exp((3p) and summed over N. The ensemble average of 
an operator A in the grand canonical ensemble is 
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3.2 Bose Systems 



3.2.1 Ideal Bose Gas 

Now we shall consider the properties of an ideal gas of particles with integer spin 
and mass m > 0. As usually for ideal gases, we neglect the interaction between 
the particles. The Hamiltonian of the ideal Bose gas is: 



H 



^ 2m v 

v 



np is the number of particles in the state with momentum p. Because we are deal- 
ing with Bosons, can also be > 1. We want to study condensation where two 
phases come into contact, the number of particles in one phase is not fixed, there- 
fore we must work with the grand-canonical distribution. The grand-canonical 
partition function of ideal Bose gas is: 



(Here we will disregard the (2S + 1) factor which comes from the spin degener- 
acy.) In the grand partition function, each sum over is over all non-negative 
integers and, since the particles do not interact, we can split the partition function 
into a product: 



s= E E 

no=0 ni=0 



no 



-/3(pl/2m-a) 



(3(p 2 2 /2m-fM) 



"2 



■}= (3.18) 




-/3(p 2 /2m- fj,) 



n- 



ze -/3p2/2m ' 



Here we have introduced the fugacity z = e^. The equation of state of an ideal 
Bose gas is (see Table I): 

PV = k B T\nE(T,V,fi) 



PV 



-J> [l-C 



-/3(p 2 /2m-M) 



(3.19) 
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The thermal average of N is determined by the chemical potential // (see Tables): 
We can write (JV) = J2pnp where 

= e /3(p2/2m- M ) _ X ( 3 - 21) 

is the Bose-Einstein distribution function, which tells the occupation probability 
of a given (non-degenerate) state p. For (N) to be finite and positive, the second 
term in the denominator has to be < 1 for any p, in particular for p = 0. This 
means that /j, must be < and < z < 1. 

For large V we replace the sum by an integral: 



and we get: 



47rl/ /*°° 2 , 1 

= ^3/2^) (™^§f) = J #3/2(2), (3-23) 



where 



1 /-oo r ^-l J 

= ^_ / dx— — = Y - (3.24) 

y K ' r(v) Jo ^z- 1 - 1 f-r Z" 



i» y 

and A is the thermal wavelength, 

/ 27Tft 2 \ V2 

A 



rnksT 



The functions # 3 / 2 (2) and #5/2(2) are shown in Fig 3.1. For later use: #3/2(1) 
2.612 and #5/2(1) 1.342. The pressure is: 

' /,; ^ (A; B T) 5 / 3 # 5 / 2 (^) = -^#5/2(2). (3.25) 



V2tt^V v " 7 y&//v 7 A3 
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Figure 3.1: The functions g 3 / 2 (z) and g 5 / 2 (z). 



Now let us start increasing the pressure at fixed (high enough) T. This increases 
the chemical potential, the fugacity, and the density of particles, (JV) /V. How- 
ever, we know that z must be < 1. What happens when z — > 1? Can we further 
increase the pressure or the density of particles? 

To answer this question we must go back to Eq. (3.20). The term with p = 
which was not included in (3.23) (The density of states n(E) has a \[E depen- 
dence on E and vanishes as E — > 0.) becomes singular (divergent) when p — > 
and it has to be treated separately. This divergence has very important conse- 
quences for a Bose gas, as we shall show now. Instead of Eq. (3.20) we must 
write: 



<ao = E 



1 



m e 



f3(p 2 /2m-n) 



V 



(3.26) 



and the equation of state is 

PV AttV 



k B T 



h 3 



p 2 dp In 



1 - ZQ -^/2m 



-ln(l -z) 




(3.27) 
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The last terms in Eqs. (3.26) and (3.27) come from the termp = and correspond 
to the particles in the lowest energy, in the ground state. In the thermodynamic 
limit the last term in (3.27) vanishes whereas it remains finite in (3.26). We will 
denote 



(N) = (N') + (N ), (3.28) 

where (N 1 ) is the average number of bosons in the continuum (excited, p > 0) 
states, 

(N') = y 3 g 3/2 (z) (3.29) 

and {N ) the average number of bosons in the ground state, 

(No) = (3.30) 

The average density of bosons in the continuum states reaches its maximum value 
when gs/2(z) is maximal, that is for z — > 1 Qx — > 0), 

^%^= 93 , 2 (l)(ff) 3/! «T^. (3.31) 

The fact that (N') is limited and cannot exceed (N' max ) has its origin in the 
quantum-mechanical and statistical nature of the system. The number of avail- 
able states in a box with volume V is finite and is equal to j^. Each state, on the 
average, accommodates n$ particles, where is given by the Bose-Einstein dis- 
tribution function, Eq. (3.21). If there are more than (N 1 ) particles in the system, 
they are pushed into the ground state. 



3.2.2 Bose Einstein Condensation 

We will now investigate the behaviour of ideal Bose gases, described by Eqs. 
(3.26, 3.27). Let us consider a system of N bosons in a volume V. 

At high enough temperature, (N' max ) is so large that iV < (N' max ) and the 
chemical potential is determined from 

N =^= V ( ! t§-f' 1 9sdz). (3.32) 
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This means that z < 1 and // < 0. The density of particles in the ground state is 

(N )/V = */(l -z)/V -> as \/^oo. (3.33) 

Almost all the particles are in the |pj > states and the ground state is macro- 
scopically empty. The equation of state simplifies to 

P _*/>(*) (334) 



k B T A 3 

At very high T, g 3 / 2 (z) and #5/2(2) are 1 and we approximate [see the series 
for n ,Eq. (3.24)]: 

#5/2(2) « 5-3/2 « 2 (3.35) 
and the equation of state becomes: 

P .Jmk B T\ m (mk B T\ m z 

z is determined by the density of particles: 

N (mk B T\ 3/2 (mk B T\ 3/2 z 

The fugacity 2 = (N)\ 3 /V vanishes as T -3 / 2 at high T and for fixed density. 
This justifies the approximation (3.35). Eliminating z from the last two equations 
yields the equation of state of the ideal Bose gas at high T: 



PV =< N > k R T 



(3.38) 



which is identical to the equation of state of a classical ideal gas (as it should be 
at high T!). 

Now we lower the temperature by keeping N and V fixed. Eq. (3.32) tells 
us that 53/2(2) must increase. This means that z and fi also increase. Eventually, 
a temperature is reached, where z reaches its maximum value (z — > 1, N = 
(N' max )). This is the transition temperature T c and is determined by the condition 

N f mk B T c \ 3/2 

v = */ 2(1 H~2^H (339) 
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from which we find: 



2nh 2 



k B m[g 3/2 (l)}^ \ V 



A\ 2 / 3 



(3.40) 



If we insert this expression into the equation of state, Eq. (3.34), we obtain the 
pressure at which the transition takes place: 



P c = k B T c /2 g 5/2 (l) 



mk 



B 



3/2 



oc T, 



5/2 



C 



(3.41) 



At Tc, the ensemble average of the density of particles in the continuum states 
reaches its maximum value whereas the ground state is still empty, (N )/V = 0. 

Below T c , (N' max ) further decreases and becomes < A (we keep A fixed). 
The chemical potential further approaches zero (/x cannot be exactly = 0, because 
that would mean that A = oo) and z — > 1 in such a way that the total number of 
particles is kept constant, z is very close to I, we can put g n (z) as g n (l) and the 
number of particles in the continuum states is 



(N 1 ) = V 



' inksT 



3/2 



^3/2(1)- 



(3.42) 



When (N') < N, the ground state starts to fill, (jV ) increases. Below T c , thus, 
only a part of the particles can be accommodated in the continuum states. The rest 
must go into the ground state! The number of bosons in the ground state is: 



(N ) — N — (N 1 ) = N - V 



rnksT 



3/2 



53/2(1) 



from which we get: 



N 



Tc) 



3/2 



(3.43) 



(A^o) /A is finite, the state with p = is occupied with a macroscopic number of 
particles, the particles condense in the momentum space into the zero-momentum 



63 



I. Vilfan 



Statistical Mechanics 



Chapter 3 




Figure 3.2: Isotherms of the ideal Bose gas. 




0.5 1 1.5 

T/Tc 



Figure 3.3: Temperature dependence of the relative number of bosons in the 
ground ((N ) /N) and excited ((TV') /N) states. 

state. This is the Bose-Einstein condensation. The fugacity is, from Eq. (3.30): 
z = (N ) /((N ) + 1). In the thermodynamic limit, when V and (N ) — > oo, 
z = 1 and /i = below T c . The equation of state is: 

h Jmk B T\ 3/ \ g 5/2 (l) (N f ) 

Only the particles in the continuum states ( in the gas phase ) contribute to pres- 
sure. The particles in the ground state (condensate) are at rest, they cannot exert 
any pressure. If the density is increased, the extra particles fall into the ground 
state and the pressure does not increase. At fixed T <T C the density of particles 
in the gas phase is constant and P is independent of V, see Fig. 3.2. The temper- 
ature dependence of (No) is shown in Fig. 3.3. In the limit T — > 0, (N' max ) — > 
and all the particles are in the ground state, N = (No). 
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p 




Gas phase 







T 



Figure 3.4: Phase diagram of an ideal Bose gas. 



The corresponding phase diagram is shown in Fig. 3.4. The normal, gas phase 
exists for T > T c , i.e., to the right of the transition line. On the line, given by 
(3.41), condensation takes place and since the particles in the ground state don't 
contribute to pressure, the condensate lies on the transition line itself. Notice that 
the "critical points" in Bose systems form a line in the P vs. T or in the P vs. V 
planes and are not in a single point like in van der Waals gases. 

We now invert Eq. (3.44), 



As T — > 0, also (N') — > and therefore V — > at any finite pressure. The 
condensed phase does not occupy any volume. This means that the ideal Bose gas 
can be compressed to zero volume without any increase in pressure. 

Some comments: The condensate (ground state) contributes neither to E, nor 
to Cy, P, or VI At low T, C v vanishes like T 3 / 2 in contrast to the photon gas, 
where Cy oc T 3 . They differ because the density of states is different in both 
cases. There are more excited states available for particles than for photons and 
the specific heat is greater. 

Of course, this, and infinite compressibility are the artifacts of the ideal gas 
model. In reality, because of atomic repulsion, the volume of the condensed phase 
does not vanish and the compressibility does not diverge. The equation of state of 
a non-ideal Bose gas is shown in Fig. 3.5. 



V = 



05/2(1) (N') 
93/2(1) P 



k B T. 



(3.45) 
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Figure 3.5: An isotherm of a non-ideal Bose gas. 




point 



T 



Figure 3.6: Phase diagram of 4 He. He I is the normal, and He II the superfluid 
phase. 



C0 A - 




Figure 3.7: Dispersion curve of elementary excitations in liquid 4 He. The excita- 
tions around the minimum at p are called rotons. 
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3.3 Liquid 4 He 



A candidate for the Bose-Einstein condensation at low T is liquid 4 He. The phase 
diagram is shown in Fig. 3.6. He I is the normal liquid and He II exhibits super- 
fluid behaviour. From the fact that the line separating the solid and the fluid phases 
is horizontal at T — > 0, we conclude with the help of the Clausius-Clapeyron 
equation, 

dp „ •S'fluid — SUiki ,~ At ~. 

^r =0= t7 t7 — (346) 

al ^fluid — ^solid 

that liquid He II at T — > has no more entropy than solid helium, which is highly 
ordered, of course. We conclude that the particles in the ground state (T — > 0) 
have already the highest possible degree of order although they are in the fluid 
and not in the solid phase. We interpret the He II phase as a mixture of two fluids, 
a normal fluid with final viscosity and a super-fluid (formed of particles in the 
ground state) with zero viscosity. Characteristic for the superfluid state is also a 
very large heat conductivity (as a consequence, He II does not boil, the atoms only 
evaporate from the surface). What is the ground state of He II and what are the 
excitations? 

We shall first discuss possible excitations above the ground state of liquid 4 He 
and only then the properties of the ground state. In the liquid, the atoms are 
close together, they are strongly interacting (exchanging momenta). A state of 
an atom with fixed momentum is not an eigenstate of the Hamiltonian, the atoms 
are subject to a collective motion. The elementary excitations in the system are 
collective excitations, quasiparticles, phonons (longitudinal because it is a liquid) 
and rotons (see Fig. 3.7). At low momenta, the excitations are phonons with a 
linear dispersion relation, 

e p = cp. (3.47) 

At high momenta (short wavelengths), the interaction between the atoms causes 
a minimum in the dispersion curve. Elementary excitations around this minimum 
are (for historical reasons) called rotons. Their energy is written in the form: 



(P ~ Po) 2 
2m* 



e p = A+ yF „_ F y , (3.48) 



where A is the energy gap, A/k B = 8.6K, and m* = 1.0 x lCT 27 kg is the 
effective mass of rotons. po = 2.0 x 10 _24 kg m/s. At low T, the density of these 
quasiparticles is low and they don't interact. 
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The grand canonical partition function for phonons (p = 0) is, in analogy 
with (3.18): 

Ins,* = J o p 2 dpln(l-e-^) 

AttV r°° x 3 dx 



3(h(3c) 3 Jo e x - 1 

- |S> T > 3 (3 - 49 > 

The grand canonical partition function for rotons (also p = 0) is: 

47rV r 00 



lnS rot - - — 



/ 00 j9 2 ^ln(l-e-' 3A - /3 ( p - p «) 2 / 2m *) 



At low T (large /?), e ^ A ^ p po) 2 / 2 " 1 * is small and we approximate ln(l — x) 
—x. (This is equivalent to using the Maxwell-Boltzmann distribution.) 



l nErot = ^ e ^ A r P 2 dpe-^ p - p0)2/2m * 
h 6 Jo 

/oo 



— e^ A /°° v 2 dve- p{ *- po)2/2 "> 



= "W~ ( P » + Tj e (3 ' 50) 

From these equations, one can calculate the specific heat (and other thermody- 
namic quantities) at low temperatures. 



3.3.1 Superfluidity 

Consider an isolated system (in which energy and momentum are conserved) com- 
posed of a long cylindrical capillary filled with liquid 4 He at T = 0. Initially, let 
the capillary tube rotate along the symmetry axis with the tangential velocity v 
and let the fluid be at rest. We would expect that - because of friction - energy 
would be transferred from the tube to the fluid and that the fluid will gradually be 
excited out of the ground state. Let us check under which conditions the energy 
and momentum can be transferred to the fluid. Immediately after the tube started 
to rotate, the tube had the kinetic energy (M is the mass of the tube) and 

the fluid was at rest and had the ground-state energy E . At a later time, we expect 
the tube velocity to decrease by Av (because we consider an isolated system) and 
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that in the fluid an (at least one) elementary excitation with energy e p and momen- 
tum p is created. The conservation laws demand that the energy at a later time 
is: 

-M(v - Av) 2 + e p + E = -Mvl + E o 

(3.51) 

and that the momentum is: 

M(v - Av) +p t = Mv =>• MAv = p t (3.52) 

where p t is the excitation momentum in the tangential direction. For Av ^ 0, the 
first equation gives the inequality 

e p = Mv Av - \m{Av) 2 < Mv Av. (3.53) 

After eliminating MAv, the conservation laws tell us that 

v-p>e p . (3.54) 

Inspection of Fig. 3.7 shows that this is only possible if v is larger than a threshold 
velocity v c . At small relative velocity between the tube and the fluid, the rotons 
are not excited, the momentum is not transferred to the fluid, the fluid will stay at 
rest. This means that there is no friction, no viscosity - this is the superfluid phase. 
(As we will see later, roton excitations are not responsible for the excitation out 
of the ground state. Vortices have higher energy, but lower e/p ratio.) 



3.3.2 Vortex Excitations 

There are other collective excitations in liquid 4 He that have higher energies than 
phonons or rotons, but which nonetheless have much lower energy to momentum 
ratio than rotons. The creation of these excitations allows energy and momentum 
to be transferred to the fluid at a much lower critical velocity than in the case of 
rotons. 

Let us now discuss such excitations of the condensate. The wave function of 
an excited collective state is 

* = ^ (r) , (3.55) 

the (mass) current density is 

-> ih 



hi/j 2 V(j)(r). (3.56) 
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On the other hand, 



A 7 



(3.57) 



3 = 



V 



mv 



where N/V = p = ip 2 . From the two equations we get for the particle velocity 

i?=— V0(r). (3.58) 
m 

The wavefunction must be single-valued and the phase is determined modulo 
2%, therefore the closed-path integral of the momentum is 



where h is the Planck's constant and n an arbitrary integer. 

To understand the nature of these excited states it is best to consider liquid 4 He 
within a small rotating bucket. If the bucket is slowly brought into rotation, the 
liquid will behave in the following way. At small angular velocities of the bucket, 
the liquid will remain at rest. If the angular velocity is increased beyond some 
threshold angular velocity, the energy is transferred to the fluid, and the fluid will 
make an abrupt transition to a state in which there is a vortex line running from 
top to bottom, parallel to the axis of rotation. The vortex line looks like a small 
whirlpool, the fluid flows around the line with the velocity that decreases with 
distance from the line. The vorticity, defined as the line integral of the momentum 
along any path encircling the vortex line, Eq. (3.59), is quantized. As the angular 
velocity is increased, more and more quantized vortex lines are created. Once a 
vortex line is created, the macroscopically occupied single-particle ground state 
is no longer a zero-momentum eigenstate. It is a quantum state with a flow pat- 
tern that includes the vortex lines. When there are many vortex lines, they will 
distribute themselves in a pattern so that the average fluid velocity (averaged over 
many vortices) is the same as in the case of rigid-body rotation. 

To summarize, the quantization of vortices is related to the phase factor of the 
condensate wave function, defined in Eq. (3.60). In the lowest-energy state, ip and 
do not depend on f, they are constant. The state in which varies with r is an 
excited state which describes the fluid with one or several vortices. 

3.3.3 Order Parameter 

A superfluid can be considered as a mixture of two coexisting fluids, of a normal 
component and of a superfluid component. At a first glance, one would choose 




(3.59) 
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the density of particles in the ground (superfluid) state as the order parameter, 
because it is 7^ in the low-temperature phase and vanishes at Tq. However, in an 
interacting system, the operator N doesn't commute with H, the single-particle 
momenta are not good quantum numbers. On the other hand we already know 

— * — * 

that the ground state is a state with zero momentum p = = — ihV V = 0. 
This means that the wavefunction is a constant in the ground state. Therefore we 
choose for the order parameter the statistical average of the ground-state wave 
function which commutes with H and which can be written in the form: 

(tf ) = - Tr ^jQ-m-^ = ^ (3 _ 60) 

(ip is the amplitude and the phase of the wavefunction). The density of the 
superfluid component is (N )/V = i\) 2 = constant. In the superfluid phase, a 
macroscopically occupied single -particle quantum-mechanical ground state with 
zero momentum extends through the whole volume available to the fluid. Each 
atom is not at any particular place, it is simultaneously everywhere in the macro- 
scopic volume! This is the meaning of the order parameter defined in (3.60). 
Symmetry breaking. There is a strong parallel between (\&) and the sponta- 

— * 

neous magnetization (M) in ferromagnets, 

(M) = Tr (Me _/W ) , (3.61) 

Zj ^ ' 

— * 

where M is the (operator for the) total magnetic moment and 

H = —J Si ■ Sj. 

In the absence of an external magnetic field, thus, H is invariant under rotations, 
M and — M occur with equal probability and therefore (M) is always zero. But 
we know that the ferromagnets do have spontaneous magnetization at low T. The 
resolution lies in "spontaneous symmetry breaking" - the ground state of a Hamil- 
tonian does not possess the full symmetry of the Hamiltonian. The ground state 
wavefunction is not rotationally invariant, because (M) points along a definite di- 
rection in space whereas H is rotationally invariant. In the case of ferromagnets, 
thus, the transition corresponds to spontaneous breaking of the rotational symme- 
try. One must therefore redefine (M) in such a way that +M and — M are not 
both included in the partition function with equal weights. This is most easily 
done by imposing an infinitesimal external field H: 

(M)= lim lim -TrMe - ^-^. (3.62) 

H^O V-+00 Z 
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Returning to Bose systems, the Hamiltonian has a global gauge invariance, 
the wave functions are invariant with respect to the phase shift, 



(3.63) 



where is an arbitrary phase factor (real number). Bose condensation corresponds 
to a spontaneous breaking of global gauge invariance. When calculating the order 
parameter, we imagine, like in magnetic systems that the system is subjected to an 
(infinitesimal) external field rj which couples linearly to the order parameter. We 
calculate the statistical average of ^(f) in the thermodynamic limit: 



(*(r)> 



lim lim 



1 



Tr*(r)e' 



-0£ 



(3.64) 



where 



£ = H- fiN - I d 3 r[ty(r)ri(f) + &(r)rjt(r)]. (3.65) 

The only essential difference with the ferromagnetic case is that the external field 
here is fictitious, it cannot be realized experimentally. 
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